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Abstract 

Superstatistics is a superposition of two different statistics relevant 
for driven nonequilibrium systems with a stationary state and inten- 
sive parameter fluctuations. It contains Tsallis statistics as a special 
case. After briefly summarizing some of the theoretical aspects, we 
describe recent applications of this concept to three different physi- 
cal problems, namely a) fully developed hydrodynamic turbulence b) 
pattern formation in thermal convection states and c) the statistics of 
cosmic rays. 



1 Introduction 



The formalism of nonextensive statistical mechanics ^1 |2l El El can be re- 
garded as an embedding of ordinary statistical mechanics into a more general 
framework. The basic idea is that some systems to which ordinary statistical 
mechanics is not applicable may still obey a similar formalism if the Shannon 
entropy is replaced by more general entropy measures. Good candidates for 
these more general entropy measures are the Tsallis entropies, which depend 
on a real parameter q and which reduce to the Shannon entropy for g 1. 
For general g, one obtains generalizations of the canonical distributions, rep- 
resented by g-exponentials, often called Tsallis distributions. There can be 
a variety of reasons why ordinary statistical mechanics is not applicable to 
a particular system: There may be long-range interactions, metastability, 
driving forces that keep the system out of equilibrium, etc. 

Tsallis distributions are indeed observed in a large variety of physical 
systems |SllllElinilZllHlini, niany of them being in a driven stationary state far 
from equilibrium. An important question is why these or similar distributions 
are that often observed in experiments. Can we give dynamical reasons 
for the occurence of Tsallis statistics in suitable classes of nonequilibrium 
systems? 

This is indeed possible. One can easily construct classes of stochastic 
differential equations with fluctuating parameters where one can rigorously 
prove that they generate Tsallis statistics (TUlilll. Indeed, for many systems 
the reason why Tsallis distributions are observed can be easily related to the 
fact that there are spatio-temporal fluctuations of an intensive parameter 
(e.g. the inverse temperature, or a friction constant, or the amplitude of the 
Gaussian white noise, or the energy dissipation in turbulent flows). If these 
fluctuations evolve on a long time scale and are distributed according to a 
particular distribution, the x^-distribution, one ends up with Tsallis statistics 
in a natural way. For other distributions of the intensive parameter, one ends 
up with more general statistics: So-called superstatistics [12,, which contain 
Tsallis statistics as a special case. Generalized entropies (analogues of the 
Tsallis entropies) can be defined for these superstatistics as well 1121 lEl , and 
generalized versions of statistical mechanics can be also constructed, at least 
in principle. It has been shown that the corresponding generalized entropies 
are stable |T3 UHl • 

In this paper we will briefly review the superstatistics concept, and then 
show that the corresponding stationary probability distributions are not only 
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a theoretical construct but of practical physical relevance. We will concen- 
trate onto three physically relevant examples of applications: Hydrodynamic 
fully developed turbulence, defect motion in convection states, and the statis- 
tics of cosmic rays. 

2 What is superstatistics? 
2.1 The basic idea 

Let us give a short introduction to the 'superstatistics' concept Con- 
sider a driven nonequilibrium systems with spatio-temporal fluctuations of 
an intensive parameter f3. This may be the inverse temperature, or a chem- 
ical potential, or a function of the fluctuating energy dissipation in the flow 
(for the turbulence application). Locally, i.e. in spatial regions (cells) where 
P is approximately constant, the system is described by ordinary statisti- 
cal mechanics, i.e. ordinary Boltzmann factors e~^^, where E is an effective 
energy in each cell. In the long-term run, the system is described by a 
spatio-temporal average over the fluctuating p. In this way one obtains a 
superposition of two statistics (that of P and that of e~^^), hence the name 
'superstatistics'. One may define an averaged Boltzmann factor B{E) as 

POO 

B{E) = / /(/3)e-^^ci/3, (1) 

where /(/5) is the probability distribution of (3. For so-called type-A su- 
perstatistics, one normalizes this effective Boltzmann factor and obtains the 
stationary long-term probability distribution 

PiE) = ^B{E), (2) 

where 

Z= B{E)dE. (3) 
Jo 

For type-B superstatistics, one includes the /^-dependent normalization con- 
stant into the averaging process. In this case 

POO 1 
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where Z{f3) is the normahzation constant of e~^^ for a given /3. Both ap- 
proaches can be easily mapped into each other, by defining a new probabihty 
density = c-/(/3)/Z(/5), where c is a normahzation constant. It is obvi- 
ous that Type-B superstatistics with / is equivalent to type-A superstatistics 
with /. 

As will be proved later, p{E) becomes a Tsallis distribution if /(/?) is 
chosen as a x^-distribution. But at the moment we keep general. Many 
relevant concepts of superstatistics can be formulated for general /(/?). 



2.2 Dynamical realization 

A simple dynamical realization of a superstatistics can be constructed by con- 
sidering stochastic differential equations with spatio-temporally fluctuating 
parameters JT] . Consider a Langevin equation for a variable u 

u = -fF{u) + crL{t), (5) 

where L{t) is Gaussian white noise, 7 > is a friction constant, a describes 
the strength of the noise, and F(u) = —-^V{u) is a drift force. If 7 and 
a are constant then the stationary probability density of u is proportional 
to e"^^*^"-*, where P ■= can be identified with the inverse temperature 
of ordinary statistical mechanics. Most generally, however, we may let the 
parameters 7 and a fluctuate so that (3 = ^ has probability density 
These fluctuations are assumed to be on a long time scale so that the system 
can temporarily reach local equilibrium. In this case one obtains for the 
conditional probability p{u\(3) (i.e. the probability of u given some value of 

P(«l/5) = ^exp{-/?V(«)}, (6) 

for the joint probability (3) (i.e. the probability to observe both a certain 
value of u and a certain value of (3) 

p{u,(3)=p{u\f3)f{(3) (7) 

and for the marginal probability p{u) (i.e. the probability to observe a certain 
value of u no matter what f3 is) 

POO 

p{u) = / p{u\p)fmp. (8) 

Jo 
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This marginal distribution is the generahzed canonical distribution of the 
superstatistics considered. The above formulation corresponds to type-B 
superstatistics. 

Let us now consider a few examples of possible superstatistics, by consid- 
ering different examples of distributions Note that f3 lives on a positive 
support, so Gaussian distributions of (3 are unsuitable. 



2.3 x^"Superstatistics 

One of the most natural choices of the probabihty density of /? is given by 
the x^-distribution 

n 

Note that if one adds Gaussian random variables then the sum is again 
a Gaussian random variable. But if one adds Gaussian random variables 
squared, one gets a x^-distributed random variable. In this sense the x^- 
distribution (also called F distribution) is a typical distribution for positive 
random variables that naturally arises in many circumstances. Let us denote 
the n independent Gaussian random variables as Xj, i = 1, . . . ,n and assume 
they have average 0. Then 

n 

P:=J2xf (10) 

i=l 

has the probability density ©. The average of the fluctuating P is given by 

/>oo 

{/3) = n{X') = / (3fiP)d/3 = (3o (H) 
Jo 

and the variance by 

{p')-p'o = IpI (12) 

For linear drift forces F{u) = —u one obtains for the conditional probability 

p{u\P) _ 



P(^l/5) = Vf^exp<^--/3^i^[>, (13) 
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and for the marginal probability p{u) a short calculation yields 

Piu) = r piu\/3)fi/3)d/3 = ^^^^7^ (^) ' T^—^- (14) 

Thus the stochastic differential equation (0) with x^-distributed P = 
generates the generalized canonical distributions of nonextensive statistical 
mechanics ^ 

p{u) ~ ^ (15) 

(1 + i/3(g - IW) 

provided the following identifications are made: 

1 1 2 

- + - ^ q = l + (16) 



g- 1 2 2 n+1 

kq-l)P=- ^ P = ^f3o. (17) 
2 n 3 — g 

2.4 Log-normal superstatistics 

Of particular interest for models of hydrodynamic turbulence is the log- 
normal distribution 

It yields yet another possible superstatistics (see [T71 UHl HH 1201 for related 
turbulence models); m and s are parameters. The average (3q of the above 
log-normal distribution is given by Po = rriy/w and the variance by cr^ = 
m?w{w — 1), where w := . One obtains for linear drift forces F[u) = —u 
the superstatistics distribution 

p{u) = — /""rf/5/5-V2expizi!^lll!u-^^"'. (19) 



27rs Jo 1 

The integral cannot be evaluated in closed form, but it can be easily numer- 
ically evaluated. 
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2.5 F-superstatistics 



Another example is superstatistics based on F-distributions. The F-distribution 
is given by 

ff3) = n{v+n^)/2) fkv_\ ^'^ pj:^ 

T{v/2)T{w/2)\w ) (l + ^/3)(-+-)/2- ^ ^ 

Here w and v are positive integers and 6 > is a parameter. We note that for 
V = 2 we obtain a TsaUis distribution. However, this is a Tsalhs distribution 
in /?-space, not in _E-space. 
The average of (3 is given by 

= (21) 

b[w — 2) 



and the variance by 



2 2w'^{y + w - 2) 



b^v{w-2Y{w~Ay ^^^^ 

Again, the integral leading to the marginal distribution p{u) cannot be 
obtained in closed form, but is easily numerically evaluated. Superstatistics 
based on F-distributions has been studied in more detail in j^I], and possible 
applications in plasma physics were sketched there. 



2.6 General properties of superstatistics 

For small all superstatistics have been shown to have the same first- 
order corrections to the Boltzmann factor of ordinary statistical mechanics 
as Tsallis statistics has [12]. For moderately large one often observes 
similar behaviour as for Tsallis statistics (see [20] for examples). On the 
other hand, the extreme tails of p{E) for very large E are very different for 
the various superstatistics. Tsallis distributions decay with a power law for 
E ^ oo, general superstatistics can have all kinds of asymptotic decays. 

For any superstatistics, not only Tsallis statistics, one can generally define 
a parameter q by the relation 

(g - l)/52 = a', (23) 

or equivalently 



For x^-superstatistics, this q coincides with the index q of the Tsalhs en- 
tropies. For the log-normal distribution we obtain from eq. (j24j) q = w and 
for the F-distribution q = 1 + The meaning of this generally de- 

fined parameter q is that ^/q'^^ = is just the coefficient of variation of 
the distribution defined by the ratio of standard deviation and mean. 

If there are no fiuctations of f3 at all, we obtain q = 1, i.e. ordinary statistical 
mechanics. Our formula (j24j) relating q and the variance of the f3 fluctuations 
is valid for both type- A and type-B superstatistics, just that one has to form 
averages with either / and /. 

As recently shown by Tsallis and Souza Olini, general superstatistics 
maximize more general classes of entropy-like functions subject to suitable 
constraints. This can be used as a starting point to construct a generalized 
statistical mechanics for general superstatistics. 

2.7 Many particles 

So far our dynamical realization of a superstatistics in terms of a stochastic 
differential equation was written down for one test particle in one dimension. 
To generalize to particles in d space dimensions, we may consider coupled 
systems of equations with fluctuating friction forces, as given by 

Ui = -'yiFi{ui, . . . ,un) + (TiLi{t) i = l,...,N. (25) 

Suppose that a potential V{ui, . . . ,itN) exists for this problem such that 
Fi = —-j^V. We can then proceed to marginal stationary distributions in a 
similar way as before. One just has to specify the statistics of the in the 
vicinity of each particle i. 

If all Pi = \ are given by the same fluctuating random variable A = /3 

i 

the integration yields marginal distributions of the form 

p(ni, . . . , n^) = ^^^-mu.,...,u^)dp (26) 

(type-B superstatistics). The result depends on both, /(/?) and Z{P). 

Let us consider an example. We consider a partition function Z{f3) whose 
/3-dependence is of the form 

Z{P) = [ dui--- duNC-^^ ~ /5^e-^^ (27) 
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where x and y are some suitable numbers. We consider x^-superstatistics for 
this problem. The marginal distributions are obtained as 



where 



and 



. . . , vTn) ~ z — (28) 



(29) 



n — 2x 



Is the assumption of a single fluctuating (3 realistic for many particles? 
No. In many physical applications, the various particles will be dilute and 
only weakly interacting. Hence in this case (3 is expected to fluctuate spa- 
tially in such a way that the local inverse temperature (3i surrounding one 
particle i is almost independent from the local (3j surrounding another parti- 
cle j. Moreover, the potential is approximately just a sum of single-particle 
potentials V{ui, . . . , un) = X^ili ^s{ui). In this case integration over all /3i 
leads to marginal densities of the form 

P(«i, ...,un)= Tdf],--- r d(3^ f[ JM^e-^^y^^'^^ (31) 

^0 ^0 ^s\Pi) 



and for x^-superstatistics one obtains 



N 



p(Mi,...,Mjv) ~ TT z ^- —. (32) 

t\{l + (i{q-l)Vs{ui)) — 

This means the A^-particle nonextensive system reduces to products of 1- 
particle nonextensive systems. 

We notice that for many particles there is no unique answer to what the 
physically relevant generalized canonical distributions are — it depends on 
the spatio-temporal statistics of the Pi. Systems of very high particle density 
may be better described by (|^. dilute systems better by (jHH). Our physical 
applications in section 3 mainly relate to the form (jHT|) . 
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2.8 Further generalizations 



We just briefly mention some further generalizations, without working them 
out in much detail. 

• We may regard the temperature T rather than p — as the funda- 
mental variable, and can then take for T a x^-distribution, log-normal 
distribution, F-distribution, etc. This case is easily reduced to the 
framework studied so far, since any chosen probability distribution 

of T implies a corresponding probability distribution of (3 given by 

fM-fT{T)\^\-Mp-')j,. (33) 

• For each particle there may be several different intensive parameters 
that fluctuate, for example not only the inverse temperature f3 but also 
the chemical potential may fluctuate in a spatio-temporal way. In 
order to proceed to the marginal distribution, one then simply has to 
do two integrals (over /3 and /i) rather than one. 

• Rather than starting from ordinary Boltzmann factors e^^^ and averag- 
ing these over a fluctuating j3, we could also construct a superstatistics 
where the statistics in the single cells is already given by Tsallis statis- 
tics, i.e. one would be /3-avcraging g-cxponcntials e^^^ in this case. This 
is related to the case of two fluctuating parameters discussed above: We 
may think of cases where the integration over the fluctuating // yields 
g-exponentials, which then still have to be averaged over (3. 

3 Applications 

3.1 Fully developed hydrodynamic turbulnce 

In the turbulence application, u stands for a local velocity difference in the 
turbulent flow. On a very small time scale, this velocity difference is essen- 
tially the acceleration. The basic idea is that turbulent velocity differences 
locally relax with a certain damping constant 7 and are at the same time 
driven by rapidly fluctuating chaotic force differences, which in good approx- 
imation are modelled by Gaussian white noise. One knows that in turbulent 
flows the energy dissipation fluctuates in space and time. Moreover, it is 
known since the early papers by Kolmogorov in 1962 that the probability 
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density of energy dissipation is approximately log-normal in a turbulent flow. 
Hence, if /3 is a simple power-law function of the energy dissipation, this im- 
plies a lognormally distributed (3. We thus end up in a most natural way 
with log-normal superstatistics. 

Fig. 1 shows an experimentally measured p{u) of velocity differences u at 
scale r = 92.5 r/ in a turbulent Taylor- Couette flow, as obtained by Swinney 
et al. 0. Tj denotes the Kolmogorov length scale. The data have been rescaled 
to variance 1. The dashed line is the prediction of log-normal superstatistics. 
Apparently, there is excellent agreement between the measured density and 
log- normal superstatistics as given by eq. (fT^ . The fitting parameter for this 
example is = 0.28. This is the only fitting parameter, m = y/w is fixed by 
the condition of variance 1 of p{u). 

Fig. 2 shows the measured probability density of the acceleration of a 
Lagrangian test particle in a turbulent flow as obtained in the experiment 
of Bodenschatz et al. [221 12^1 ■ Log-normal superstatistics with = 3.0 
yields a good fit. Since Bodenschatz's data reach rather large accelerations a 
(in units of the standard deviation), the measured tails of the distributions 
allow for a sensitive distinction between various superstatistics. Though 
X^-superstatistics (= Tsallis statistics) yields also a reasonably good fit of 
the data for moderately large accelerations, for the extreme tails log-normal 
superstatistics seems to do a better job. The main difference between x^- 
superstatistics and log-normal superstatistics is the fact that p{a) decays with 
a power law for the former ones, whereas it decays with a more complicated 
logarithmic law for the latter ones. 

3.2 Chaotic defect motion in inclined layer convection 

Let us now consider another physically relevant example, so-called 'defect tur- 
bulence'. Defect turbulence shares with ordinary turbulence only the name, 
otherwise it is very different. It is a phenomenon related to convection and 
has nothing to do with fully developed hydrodynamic turbulence. Consider 
a Raleigh-Benard convection experiment: A liquid is heated from below and 
cooled from above. For large enough temperature differences, interesting 
convection patterns start to evolve. An inclined layer convection experiment 
[71 121] is a kind of Raleigh-Benard experiment where the apparatus is tilted 
by an angle (say 30 degrees), moreover the liquid is confined between two 
very narrow plates. For large enough temperature differences, the convection 
rolls evolve chaotically. Of particular interest are the defects in this pattern. 



10 



i.e. points where two convection rolls merge into one (see Fig. 3). These de- 
fects behave very much like particles. They have a well-defined position and 
velocity, they are created and annihilated (in pairs of course), and one can 
even formally attribute a 'charge' to them: There are positive and negative 
defects. For more details, see [71 121] and references therein. 

The probability density of defect velocities has been quite precisely mea- 
sured. As shown in Fig. 4, it quite precisely coincides with a Tsallis distribu- 
tion with q ^ 1.46. Also other dynamical properties of the defects coincide 
quite well with typical predictions of nonextensive models jTj. In total, these 
defects seem to behave very much like an ideal gas of nonextensive statistical 
mechanics with entropic index q 1.5. Defect turbulence thus may serve as 
an example where generalized versions of statistical mechanics are not only 
mathematically beautiful but physically useful. Note that this is an appli- 
cation to a physical system far from equilibrium, where ordinary statistical 
mechanics has little to say. 

Our dynamical realization in terms of the generalized Langevin equation 
with F{u) = —u and fiuctuating effective friction 7 makes sense as a very 
simple model for the defect velocity u. While ordinary Brownian particles 
have constant damping due to Stokes' law 7 = where v is the kinematic 

viscosity of the liquid, p is its density, m is the mass of the particle and a is 
the radius of the particle, defects are no ordinary particles: they have neither 
a well-defined mass m nor a well-defined radius a and thus one expects that 
there is an ensemble of damping constants which depend on the topology of 
the defect and its fiuctuating environment. In particular, the fastest velocities 
result from circumstances in which the defect is moving in a local environment 
with only a very small effective damping 7 acting. The driving forces L{t) are 
hardly damped during such a time interval, and lead to very large velocities 
for a limited amount of time, until another region with another 7 is reached. 

The experimentally observed value q ~ 1.46. ..1.50 for the defect statistics 
means according to eq. ()16|) that there are effectively about three independent 
degrees of freedom that contribute to the fluctuating local defect environ- 
ment. We do not know where these three effective degrees of freedom come 
from, but a very simple picture would be that the fluctuating environment of 
the defect is mainly characterized by the states of the three convection rolls 
that merge when forming a defect. 
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3.3 Statistics of cosmic rays 

Our third example is from high energy physics. We will proceed to extremely 
high temperatures, where (similar as in defect turbulence) particles are cre- 
ated and annihilated in pairs. Nonextensive statistical mechanics has been 
shown to work well for reproducing experimentally measured cross sections 
in e^e~ collider experiments |8i|25j. Here we apply it to high energy collision 
processes as induced by astrophysical sources, leading to the creation of cos- 
mic ray particles j^EI that are ultimately observed on the earth. The idea to 
apply nonextensive statistics to the measured cosmic ray spectrum was first 
presented in Pj , some related work can also be found in [^Zl |2H1 • 

Experimental data of the measured cosmic ray spectrum [^H] are shown in 
Fig. 5. Also shown is a curve that corresponds to a prediction of nonextensive 
statistical mechanics. Up to energies of 10^^ eV, the measured flux rate 
of cosmic ray particles with a given energy is well fitted by a generalized 
canonical distribution of the form 

piE) = C . — . (34) 

This is an F-distribution in E. E is the energy of the particles. 



E = J c^pI + c^Py + c^pI + m^c*!, (35) 



P = T~i is an effective inverse temperature variable, and C is a constant 
representing the total flux rate. For relativistic particles the rest mass m can 
be neglected and one has E ^ c\p\. The distribution (jMjl is a g-generalized 
relativistic Maxwell-Boltzmann distribution in the formalism of nonextensive 
statistical mechanics. The Tsallis distribution is multiplied with E'^, taking 
into account the available phase space volume. As seen in Fig. 5, the cosmic 
ray spectrum is very well fitted by the distribution ()34|) if the entropic index 
is chosen as g = 1.215 and and if the effective temperature parameter is given 
by f = /3-1 = 107 MeV. 

The above effective temperature is of the same order of magnitude as the 
so-called Hagedorn temperature Th PUl , an effective temperature well known 
from collider experiments. The Hagedorn temperature is much smaller than 
the center-of-mass energy Ecus of a typical collision process and represents 
a kind of 'boiling temperature' of nuclear matter at the confinement phase 
transition. It is a kind of maximum temperature that can be reached in a 
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collision experiment. Even largest Eqms cannot produce a larger average 
temperature than due to the fact that the number of possible particle 
states grows exponentially. The Hagedorn theory of scattering processes is 
known to work well for energies Ecms < 10 GeV, whereas for larger energies 
there is experimental evidence from various collision experiments that power- 
law behaviour of differential cross sections sets in. This power-law is not 
contained in the original Hagedorn theory but can be formally incorporated 
if one considers a nonextensive extension of the Hagedorn theory 

Let us now work out the assumption that the nonextensive behaviour 
of the measured cosmic ray spectrum is due to fluctuations of temperature. 
Assume that locally, in the creation process of some cosmic ray particle, some 
value of the fluctuating inverse temperature P is given. We then expect the 
momentum of a randomly picked particle in this region to be distributed 
according to the relativistic Maxwell-Boltzmann distribution 

Pirn = zip-f'^''''- (36) 

Here p{E\f3) denotes the conditional probability of E given some value of f3. 
We neglect the rest mass m so that E = c\p\. The normalization constant is 
given by 

2 

Now assume that f3 is x^-distributed. The observed cosmic ray distribution 
at the earth does not contain any information on the local temperature at 
which the various particles were produced. Hence we have to average over all 
possible fluctuating temperatures, obtaining the measured energy spectrum 
as the marginal distribution 

POD 

p{E) = / p{E\P)fiP)dp. (38) 
Jo 

The integral (|HH|) with /(/5) given by Q and p{E\P) given by (jSHI) is easily 
evaluated and one obtains eq. with 

5 = 1 + ^ (39) 
n + b 

and 



Z(P) = / E^-^'^dE = ± (37) 
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This just corresponds to eqs. and (jHUI) for x = — 3 and y = 0. 

The variables Xj in eq. (llUj) describe the independent degrees of freedom 
contributing to the fluctuating temperature. At very large center of mass 
energies, due to the uncertainty relation r ~ 1/Ecms, the probed volume 
is very small, and all relevant degrees of freedom in this small volume are 
basically represented by the 3 spatial dimensions into which heat can flow. 
We may physically interpret Xf as the heat loss in the spatial i-direction, 
i = x,y, z, during the collision process that generates the cosmic ray particle. 
The more heat is lost, the smaller is the local T, i.e. the larger is the local /3 
given by eq. (fTUj). The 3 spatial degrees of freedom yield n = 3 or, according 

to eq. dnni), 

g=H = 1.222. (41) 

For cosmic rays Eqms is very large, hence we expect a g-value that is close 
to this asymptotic value. The fit in Fig. 5 in fact uses q = 1.215, which 
agrees with the predicted value in eq. (PT|) to about 3 digits. It also coincides 
well with the fitting value q = 1.225 used by Tsallis et al. |9j using multi- 
parameter generalizations of nonextensive canonical distributions. 

For smaller center of mass energies, the volume probed will be bigger 
and more effective degrees of freedom within this bigger interaction region can 
contribute to the fluctuating temperature. Hence we expect that for decreas- 
ing EcMS ^ will increase and q will become smaller than 11/9. Experimental 
data measured in e^e~ annihilation experiments are in good agreement with 
the following parametrization (2b\ : 

^(^CMs) = 9 ^ (42) 

with Eq ^ 45.6 GeV. Solving for Eqms we get 

11 - 9o 

^CAfs = -^olog^^, (43) 
1 + g 

and putting in the fitted value q = 1.215 of the cosmic ray spectrum we 
can give an estimate of the average center of mass energy of the ensemble of 
astrophysical accelerators: 

EcMS ~ 161 GeV. (44) 

It should be clear that this is a very rough estimate only, moreover, it is an 
average over the unknown ensemble of all accelerating astrophysical sources. 
Some astrophysical objects can definitely accelerate to larger energies. 
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The 'knee' and 'ankle' in the cosmic ray spectrum, occuring at ~ 10^^ 
eV and E ^ 10^^ eV, are produced by other effects which can be easily 
embedded into this nonextensive model (see [2H] for details). The advantage 
of our superstatistics model is that a concrete prediction for q can be given. 
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Fig. 1 Histogram of velocity differences u as measured in Swinney's exper- 
iment and the log-normal superstatistics prediction eq. (|19|) with = 0.28. 
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Fig. 2 Histogram of accelerations a as measured in Bodenschatz's exper- 
iment, the log-normal superstatistics prediction eq. (fT^ with = 3.0, and 
Tsallis statistics with q = 1.5. 
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Fig. 3 A positive defect (black box) and a negative defect (white box) as 
seen in the experiment of Daniels et al. |24j . 
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Fig. 4 Measured distribution of positive and negative defect velocities 
and comparison with a Tsallis distribution with q = 1.46. The dimensionless 
temperature difference is ^ — 1 = 0.17. All distributions have been rescaled 
to variance 1. 
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Fig. 5 Measured distribution of primary cosmic ray particles with a given 
energy E and comparison with the distribution ()34|) with q = 1.215 and 
p-^ = 107 MeV (sohd line). 
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